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We suggest that if a localized phase at nonzero temperature T > exists for strongly disordered 
and weakly interacting electrons, as recently argued, it will also occur when both disorder and 
interactions are strong and T is very high. We show that in this high-T regime the localization 
transition may be studied numerically through exact diagonalization of small systems. We obtain 
spectra for one- dimensional lattice models of interacting spinless fermions in a random potential. 
As expected, the spectral statistics of finite-size samples cross over from those of orthogonal random 
matrices in the diffusive regime at weak random potential to Poisson statistics in the localized regime 
at strong randomness. However, these data show deviations from simple one-parameter finite-size 
scaling: the apparent mobility edge "drifts" as the system's size is increased. Based on spectral 
statistics alone, we have thus been unable to make a strong numerical case for the presence of a 
many-body localized phase at nonzero T. 

PACS numbers: 



I. INTRODUCTION 

Although Anderson's original paper on localization 
is mostly remembered for its ground-breaking re- 
sults about single particles in random potentials, one 
goal of that paper was to learn about transport proper- 
ties of highly-excited many-body eigenstates, e.g. quan- 
tum diffusion of nuclear moments. This latter goal was 
mostly neglected in subsequent research on localization 
and metal-insulator transitions. However, these ques- 
tions have been recently brought to our attention by 
Basko, et al. 0, who present detailed arguments that in- 
teracting electrons in static random potentials can have a 
true metal-insulator transition at a nonzero critical tem- 
perature. Thus these systems are argued to have an 
insulating phase, with strictly zero ohmic conductivity, 
even at a nonzero temperature. For some work on these 
questions published before Basko, et at, see for example 

In practice, few transport measurements are possible 
without first equilibrating the sample with its environ- 
ment in order to establish a steady state (by removing 
Joule heat). In metals this coupling to the environment, 
provided it is not too strong, does not affect the conduc- 
tivity (non-linear transport is another story altogether, 
see e.g. |l0|). In Anderson insulators, however, the 
heat bath plays a far less subtle role: it is what permits 
transport. Conduction occurs by variable-range-hopping, 
which is an inelastic process requiring a heat bath that 
can locally supply or absorb the energy needed to permit 
hopping of the charge carriers between localized states 
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that are not precisely degenerate. At the heart of this 
extreme sensitivity of the dynamics of a localized insula- 
tor to the coupling with its environment is its inability 
to self-equilibrate. It is therefore useful to turn the issue 
around by distinguishing conductors from true T > in- 
sulators by whether the many-particle system itself con- 
stitutes a heat bath. For example, one might ask whether 
external local probes can deposit limitless amounts of en- 
ergy or if they tend to saturate the spectrum. Similarly, 
whether or not attached leads themselves can effectively 
remove heat from the sample will generally depend on 
heat conductivity of the sample itself. Thus we see that 
whether or not a quantum system of many interacting 
degrees of freedom constitutes a heat bath is not only a 
very fundamental question, but also one of some practical 
relevance. 

To the extent that one of the most successful theo- 
ries of nature, namely thermodynamics, is founded on 
the assumption of ergodicity, we expect true insulators 
(where this assumption is strongly violated) to be rare 
and require fine tuning of some sort. The noninteracting 
Anderson insulator is one example, where the unrealistic 
condition of no interparticle interactions is crucial. Re- 
markably, the authors of Ref. 2 argue that a nonzero 
temperature Anderson insulator can be stable against 
the dephasing effects of interparticle interactions, mak- 
ing this state a sufficiently realistic possibility to be taken 
seriously and looked for in experiments (provided deco- 
herence from the rest of the universe can be ignored to a 
good approximation). 

The calculations of Ref. 2 are based on a low energy 
effective Hamiltonian whose connection with the param- 
eters of the original model of interacting electrons in a 
random potential could not be established analytically. 
Thus, it is interesting and likely worthwhile to test their 
results using other methods, and to try to learn more 
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about the nature of the proposed T > diffusive-to- 
insulating phase transition and about the range of mod- 
els that may exhibit it. We report here on one such at- 
tempt. To start we observe that application of the quan- 
titative estimates of the localization transition in Ref. 2 
to a lattice model with finite entropy and energy densities 
(i.e. finite number of states at each site) implies that the 
aforementioned localized phase and, therefore, the phase 
transition to the diffusive state can persist all the way to 
infinite temperature. This seemingly innocuous obser- 
vation has at least two important practical implications. 
First, by adapting familiar high temperature expansion 
techniques we can more or less rigorously rule out the 
possibility that such a transition is accompanied by a 
thermodynamic signature both at infinite temperature 
and by continuity at any finite temperature 11]. Per- 
haps more interestingly, the very large (exponential in 
volume) number of states available to the system at high 
temperatures can sometimes create favorable conditions 
for quickly approaching the thermodynamic limit in vari- 
ous thermodynamic and dynamic quantities, which raises 
the possibility of looking for the signs of this physics nu- 
merically, e.g. in exact spectra of finite samples. Our 
choice of the model and method of analysis are summa- 
rized below, followed by results and some preliminary 
conclusions. 



II. ENSEMBLE OF HAMILTONIANS 

To reach as large as possible a distance with a given 
size many-body Hilbert space, we study spinless fermions 
hopping and interacting on a one-dimensional lattice of 
L sites with a random potential and periodic boundary 
conditions. This model has only two states (empty and 
occupied) per lattice site. The Hamiltonian is 

H = ^K«i + V{n t - ~)(ni+i - i) 

i 

+c\c l+1 + c\ +1 c t + c\c l+2 + 4 +2 Ci\ . (1) 

The nearest-neighbor interaction is chosen to be V = 2, 
although we have explored other values. The hopping 
matrix elements to both nearest and second-neighbor 
sites are chosen to be t = t' = 1, although again we 
have explored other values. The second-neighbor hop- 
ping is included so that the model remains nonintcgrablc 
(quantum chaotic) and thus diffusive at zero randomness 
|l2j. The on-site potentials Wi are independent Gaus- 
sian random numbers with mean zero and variance W 2 . 
Each realization of the disorder potential will generally 
have mean-square random potential Yli=i w fl-^ that is 
not precisely W 2 . We have found that restricting our 
ensemble of samples to those with mean-square random 
potential precisely W 2 reduces our statistical uncertain- 
ties by about a factor of 2 in the largest samples. This 
change of statistical ensembles does produce quantita- 
tive changes in the spectral properties (mostly notice- 



able for intermediate values of disorder, 4 < W < 9), but 
it does not appear to produce any qualitative changes 
in the finite-size scaling behavior that would affect our 
conclusions and it cannot affect the system's intensive 
properties in the thermodynamic limit. 

We study all many-body eigenenergies of this Hamilto- 
nian, weighting them equally; the data shown here are for 
half- filling, L /2 particles. Thus we are studying tempera- 
tures high compared to the energy scales of this Hamilto- 
nian. If a localized phase does indeed exist in this model, 
it should be present even at high T for strong enough 
disorder. An important motivation for this choice of a 
model was our recent work |l2l ] on the same model in 
the absence of randomness, where the approach to ther- 
modynamic limit was rapid enough to observe the onset 
of hydrodynamic behavior with < 9 particles. Here we 
are limited to somewhat smaller sizes, since the random 
potential violates momentum conservation; we focus on 
sizes up to L — 16. The number of realizations needed 
to achieve adequate statistical certainty depends strongly 
on W and even more so on L. At L = 8 we average over 
10, 000 realizations whereas at L — 14, 16 only about 50 
to 100 suffice except in the putative critical region, where 
we average over 1000 realizations for each W. 



III. METHOD OF ANALYSIS 

To look for the diffusive-to-insulating phase transition 
in this model, we have chosen to use what appears to 
be numerically the most accessible quantity that shows a 
clear, well-understood difference between the two phases, 
namely the spectral statistics of adjacent energy levels of 
the many-body Hamiltonian. In the localized, insulating 
phase (assuming it exists in our many-body system), in 
the thermodynamic limit of a large sample, the eigen- 
states are localized in the many-body Fock basis of local- 
ized single-particle orbitals, so states that are nearby in 
energy are far apart in this Fock space and do not inter- 
act or show level repulsion. As a result, nearby energy 
levels are simply Poisson distributed . In the diffusive 
phase, on the other hand, the level statistics of a large 
sample are those of random matrix theory, the Gaus- 
sian orthogonal ensemble (GOE) in particular. For the 
finite-length samples that we can diagonalize, the level 
statistics cross over smoothly between these two limit- 
ing behaviors as the strength of the random potential is 
varied. This crossover becomes sharper as the length L 
is increased, and we can look for a phase transition us- 
ing standard finite-size scaling techniques; this approach 
works well for the single-particle localization transition 
in three dimensions (see, e.g., 0]). 

The choice of a quantity to compute and use for the 
finite-size scaling analysis is to some degree arbitrary: 
the hypothesis of universality implies that many features 
of the distribution of eigenvalues of the Hamiltonian are 
universal in the thermodynamic limit |l3j . By analogy to 
the Binder ratio for phase transitions with a local order 
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parameter [15| , we seek a dimensionless measure of spec- 
tral statistical properties, say r(W,L), that is expected 
to take different finite values in the thermodynamic limit, 
L — » oo, in the two phases and at the critical point 
(W > W c , W < W c and W = W c ). Since the zero of 
energy is arbitrary, it is natural to work with gaps be- 
tween many-body levels. Here in particular we consider 
gaps between adjacent many-body levels, 

S n = E n+ i — E n > , 

where the eigenvalues of a given realization of the Hamil- 
tonian for a given total number of particles, {E n }, are 
listed in ascending order. 

The dimensionless quantity we [l6| have chosen to 
characterize the correlations between adjacent gaps in 
the spectrum is the ratio of two consecutive gaps 

< r n = min{<5„, 5 n -i}/ max{£„, 5 n -i} < 1 ■ 

For uncorrelated Poisson spectrum the probability distri- 
bution of this ratio r is Pp(r) = 2/(1 + r) 2 , and its mean 
value is (r)p = 2 In 2 — 1 = 0.386. The numerically- 
determined probability distribution [17| for large GOE 
random matrices is shown in Fig. 1; its mean value 
is (r) G OE = 0.5295 ± 0.0006. As expected, level re- 
pulsion/spectral rigidity in the GOE spectra manifests 
itself in the vanishing of the probability distribution 
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FIG. 1: (color online) Disorder averaged probability distri- 
bution, P(r), for Poisson (solid black line) and GOE dis- 
tributed eigenvalues [l7l | (black dots) and for our interact- 
ing fermion model at randomness W — 3 (green, diffusive 
regime), W — 11 (red, localized regime) and W = 7 (blue, 
intermediate) for length L — 16. 



IV. RESULTS 

With interaction and hopping terms fixed as above 
(t = t' = V/2 = 1), we vary the strength of disorder 
from W=l to TF=10 or more and for each (L, W) we 



diagonalize a large number realizations, R (see above). 
For each sample we compute the spectral average of r, 
(r), over all states. We then disorder-average this quan- 
tity, [(r)], to arrive at r(W, L) exhibited in Figure 2. The 
statistical uncertainties in [(r)] are estimated as usual as 
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FIG. 2: (color online) Size L and disorder W dependence of 
r(W,L). The curves correspond to L=8, 10, 12, 14, 16 from 
top to bottom for large W . Bottom: an enlargement of the 
crossing region to make the drift of the crossings more visible. 
Where not visible, the error bars are smaller than the points. 

As expected, larger samples have more Poisson-likc 
statistics than smaller ones for strong disorder, W > 8, in 
an apparently localized regime; while for weak disorder, 
W < 4, the level statistics converge towards GOE with 
increasing L, since this is the diffusive phase. We have 
checked that the entire probability distributions P(r) 
in these regimes approach those of Poisson and GOE 
spectra (see Fig.l). There is an additional crossover at 
very weak disorder: as crystal momentum conservation 
is recovered there appears a turnaround in the statis- 
tics as the decoupling of different momentum sectors sup- 
presses the average of r below its GOE value. This latter 
crossover at weak randomness is a nuisance for us and we 
steer clear of it the best we can by working away from 
the clean limit and also not considering very short chains 
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(with less than 8 sites) where this momentum pseudo- 
conservation persists to larger values of disorder: r(W, 8) 
shows a remnant of this crossover at W = 1 , while larger 
values of L do not show it at all over the range of W 
considered. 

The simplest one-parameter finite-size scaling scenario 
for the proposed diffusive-to-insulating phase transition 
would have these traces of r(W, L) vs. W at fixed L in 
Fig. 2 all cross at W c as L — » oo, with a slope that in- 
creases with increasing L (e.g., see Q). However, we 
find that the crossings of the r(W) curves for adjacent 
L's take place at points that, as L is increased, "drift" 
progressively towards larger W and smaller (more insu- 
lating) r; see Fig. 2. As this drift precludes the straight- 
forward quantitative analysis of our data in terms of one- 
parameter scaling theory, we have exerted considerable 
effort to attempt to eliminate it jl8f], including looking 
at other temperatures, interactions and fillings, candi- 
date scaling variables other than r, and selective spectral 
averaging (e.g. excluding states in the high and low en- 
ergy tails of the spectrum) . While this drift of the cross- 
ings can be reduced (particularly by trimming tails or 
reducing the temperature) it appears that it is intrinsic 
to this model's spectral statistics and none of the many 
things we have tried eliminated or reversed it. Accepting 
this, there are two very distinct possible implications: 
either the drift converges to a finite W c (and likely to 
r c = (r)p; see below) in the large L limit, or it continues 
indefinitely to W c — oo which would imply that the in- 
sulating phase does not exist at these high temperatures. 
In fact, this latter possibility has already been advocated 
in previous work, see e.g. ref. where it was argued 
that W C (L) ~ L (i.e. W C (L) is a length dependent scale 
at which spectral statistics changes from Poisson-like to 
GOE). 

Although at this point wc cannot choose between 
these two possibilities based on these data for the spec- 
tral statistics, it is worth making some more comments 
about the former possibility |l8|: The apparent drift of 
the crossing points, {W C (L), r c (L)}, is indeed substan- 
tial along the vertical axis, as would be expected if r c is 
converging to the Poisson-limit value (r)p. Thus these 
data seem consistent with a large L limiting behavior 
W C (L) — > W c < oo, r c (L) — > (r)p, whereby the crit- 
ical point is insulating as far as level statistics are con- 
cerned. There are independent reasons, based on analogy 
to Anderson localization on high-dimensional and Cay- 
ley graphs 0| , to expect such a behavior. The proposed 
many-body localization transition is a localization transi- 
tion in an infinite-dimensional Fock space |2|. Given that, 
there may be plenty of room in that space for the states 



at a diffusive-to-localized transition (i.e., at the mobil- 
ity edge) to have an infinite localization length but still 
have a negligible overlap between states and thus no level 
repulsion and Poisson level statistics. This would imply 
that the spectral statistics should converge to Poisson as 
L increases both within the localized phase and at the 
transition, and thus the "crossings" in our Fig. 2 must 
move down to r = {r)p in the large L limit. This sce- 
nario, with a localized phase for W > W c , seems qualita- 
tively consistent with the data we have presented above. 
Unfortunately, if this is indeed the case then spectral 
statistics are not a good tool for simple finite-size scal- 
ing analysis. We shall explore other approaches to this 
problem in the near future. 



V. SUMMARY AND OUTLOOK 

We have looked for signatures of the proposed many- 
body localization transition in the statistics of exact spec- 
tra of a one-dimensional tight-binding model of strongly- 
interacting spinless fermions in a random potential. Al- 
though some indications of this phase transition are 
clearly seen, there are rather strong deviations from 
and/or corrections to finite-size scaling present. The lat- 
ter might be interpreted as calling in to question the ex- 
istence of the proposed many-body localized phase at the 
high temperatures we study. Alternatively, this failure of 
simple one-parameter finite-size scaling might be because 
the critical point has insulator-like spectral statistics. 

In closing, it may be worth noting that thus far we have 
focussed on the most elementary aspects of many-body 
localization. These may not be necessarily the easiest to 
study experimentally. Finite-size effects in dynamical re- 
sponse functions, i.e. conductivity, appear more delicate 
but they are certainly worthwhile understanding, as data 
may already exist in regimes of interest, in materials as 
diverse as magnetic salts and disordered conductors. 
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